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We measure current- voltage characteristics as function of magnetic field and temperature in Nb 
strips of different thickness and width. The instability voltage of the flux flow state related to the 
vortex critical velocity v* is studied and compared with the Larkin-Ovchinnikov theory. Beside 
the usual power- law dependence v* w B -1 ^ 2 , in the low field range a new cross-over field, B cr i, is 
observed below which v* decreases by further lowering the external magnetic field B. We ascribe 
this unexpected cross-over to vortex channeling due to a fan-like penetration of the applied magnetic 
field as confirmed by magneto-optic imaging. The observation of B cr i becomes a direct evidence of 
a general feature in type-II superconducting films at low fields, that is a channel-like vortex motion 
induced by the inhomogeneous magnetic state caused by the relatively strong pinning. 

PACS numbers: 74.25. F-, 74.25.Wx, 74.78.-w 
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I. INTRODUCTION 



In type-II superconductors non linear effects of the flux 
flow dynamics of Abrikosov vortices driven by high bias 
currents have been deeply studied in the past in the par- 
ticular case of homogeneous state in weak pinning ma- 
terials, regardless the pinning mechanism of the vortex- 
defect interaction^ - — On the contrary, an attractive in- 
terest has been recently devoted to different supercon- 
ducting materials in which inhomogeneous states are in- 
duced by the type of pinning, intrinsic or extrinsic, caus- 
ing nonlinearities of the flux flow motion resulting in a 
preferential as well as a directional motion of vortices. 4 8 
At high driving currents the current-voltage (I-V) char- 
acteristics show an abrupt switching from the linear flux 
flow state into the normal state at a threshold voltage V* . 
This voltage jump can be accounted for different intrin- 
sic or extrinsic mechanisms such as for example Joule 
heatingAi^ An intrinsic non linear conductivity based 
on non equilibrium effects has been predicted by Larkin 
and Ovchinnikov (LO). 1,2 In this approach a sufficiently 
high vortex velocity v determines the escape of quasipar- 
ticles from the vortex core. The finite energy relaxation 
time r e of quasiparticles leads to appreciable changes of 
the quasiparticles distribution and to the shrinking of the 
vortex corei^ The vortex velocity dependence of the flux 
flow resistivity implies an instability of the moving vor- 
tex lattice at a threshold velocity v* with a sudden jump 
to the normal state. 

Within the LO approach, in a magnetic field range in 
which a spatially uniform quasiparticle distribution can 



be assumed 1 , a magnetic field independent critical veloc- 
ity v Jo is derived as function of the reduced temperature 
t = T/T c ± 
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Here D is the quasiparticle diffusion coefficient and ((x) 
is the Riemann zeta function. The LO prediction has 
been confirmed on several superconducting materials and 
multilayers above a cross-over field B cr , while below B cr , 
the power-law dependence w B~ x l 2 was found^ - — 
This well known cross-over effect has been justified to 
ensure the spatial homogeneity of the non-equilibrium 
quasi-particles distribution in low magnetic fields J^, 
In this paper we focus on a very low field region in wide 
Nb strips, in which we identify a new cross-over in the 
magnetic field dependence of the vortex critical velocity 
v*. We study this field regime by both I-V and Magneto- 
Optical-Imaging (MOI) measurements in order to eluci- 
date the role of the first flux penetration in the high ve- 
locity vortex dynamics. 

The observation of a new cross-over field, B cr \, from the 
regime where v* increases with decreasing magnetic field 
(v* ~ I? -1 / 2 ) to another regime where v* decreases with 
further decreasing the field, can be ascribed to an inho- 
mogenous state induced by the flux penetration. In fact, 
in the very low field region, B < B cr \, the flux does not 
penetrate with a smooth advancing front, but instead as 
a series of irregularly shaped protrusions, 16 causing the 
formation of an inhomogeneous magnetic field distribu- 
tion as observed by MO L 17 i 18 This kind of flux penetra- 
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tion is compatible with the presence of circular defects 
or micro-pores, as revealed by a FESEM analysis, which 
shows they are large enough to preferentially channel the 
moving vortices. 



II. EXPERIMENT 

We performed I-V measurements by the standard four 
contact technique on Nb films of different thickness in 
the range d = 30 -=- 150 nm. The samples were grown 
on Si(100) substrates by a UHV dc diode magnetron 
sputtering. The films were deposited at typical rates of 
0.28 nm/s. The Nb strips were obtained by standard 
photolithographic technique which provides four sam- 
ples of different widths (w =20, 40, 50, 100 jam) on the 
same substrate, with a length between the voltage con- 
tacts L = 2 mm. All samples are characterized by an 
RRR(= P3ook/piok) > 2. The Nb films are character- 
ized by a Ginzburg-Landau parameter n > 6, being the 
coherence length £o = 7 nm and the electron mean free 
path / = 5 nm with the typical value of the penetration 
depth Xl = 39 nm, 19 so they show type-II superconduct- 
ing properties. The sample parameters are summarized 
in the Table I. 

TABLE I. Physical parameters of the samples 
Sample d (nm) w (/im) T c (K) AT C (K) />iok (/i^-cm) 



NbA02 


30 


40 


7.6 


0.2 


24.0 


NbA04 


30 


100 


7.6 


0.2 


25.0 


NbAl 


60 


20 


8.2 


0.3 


15.5 


NbA3 


60 


50 


8.1 


0.2 


16.5 


NbA4 


60 


100 


8.2 


0.2 


16.5 


NbB5 


100 


100 


8.5 


0.1 


15.3 


NbC6 


135 


50 


8.7 


0.1 


16.5 


NbD8 


150 


50 


9.2 


0.2 


9.5 



We investigated a temperature range down to 0.5 T c and 
a broad magnetic field interval, perpendicularly applied 
to the film surface. The superconducting transition tem- 
perature T c is the value corresponding to R < 10 -4 ^, 
well below the thresholds used to compute the transi- 
tion width AT C , estimated by the criterion from 10% 
to 90%. The slightly high values of piox are related to 
the fabrication process. The typical value of the criti- 
cal current density at T = 4.2 K in zero external field is 
J c « 3 • 10 10 A/m 2 . In our experiments we have previ- 
ously performed a systematic study on the influence of 
self-heating by changing the bias current operation mode 
in order to ensure that the observed flux flow instabil- 
ity (FFI) is not triggered by Joule self-heating, as re- 
ported elsewhere i 19 i 2Q In order to minimize the unavoid- 
able heating effects at high driving currents, rectangu- 
lar pulses have been used with a current-on equal to 2.5 
ms and a current-off equal to 3.5 ms. We also compare 



our experimental I-V data with the numerical results ob- 
tained by J. Maza et al*2I Main differences are that the 
calculations have been performed for T close to T c with 
no applied magnetic field, and the numerical building 
procedure of the I-V curve is made up by a staircase 
ramp of current, which is not the bias operation mode 
we have experimentally performed. Instead, to minimize 
self-heating, we have chosen the short pulses mode, each 
pulse of length equal PW with a delay equal to PD, as de- 
scribed elsewhere* 2 ^ We have studied the influence of the 
pulse width and pulse delay on the I-V measurements, 20 
thus we established that PW = 2.5 ms is the maximum 
value and PD = 3.5 ms is the minimum value which are 
compatible with a reproducible I-V curve, so that the 
same instability point (I*,V*) is reached. By sweeping 
upward and downward the pulse amplitude, no hysteresis 
in I-V curves has been observed. From standard thermal 
analysis and from the magnetic field dependent dissipa- 
tion power P* = I* V* (I* is the current value at the in- 
stability point), we excluded that the observed threshold 
voltages V* were determined by heating mechanisms. 20 
Indeed, if it is the case of thermal runway, such P* should 
result magnetic field independent^ 
We also underline that I-V measurements have been per- 
formed following both Zero Field Cooling (ZFC) and 
Field Cooling (FC) procedures, leading to the same re- 
sults, as expected, since for each I-V curve, data have 
been recorded by driving the system up to the normal 
state with sufficiently high bias current. 
I-V curves of the sample NbAl, for different values of 
the magnetic field at T = 4.2 K, displaying the volt- 
age jumps, are shown in Fig. [TJ In particular the I-V 
curves for fields below B < 10.0 mT have been investi- 
gated with a fine variation of the external field around 
0.1 mT. From the measured I-V data, the critical vortex 
velocity v* (P) = V* / B L has been deduced for each mag- 
netic field value. 23 The magnetic field dependence of 
is presented in Fig. [2^ and [2)3, for the NbB5 and NbA02 
samples respectively, in the low field range B < 50.0 mT. 
Whereas for larger fields, it has been clearly observed the 
usual cross-over field B cr from the regime where is 
proportional to P -1 / 2 , to another regime at higher fields 
where v* is found almost constant ] 19 i 24 . From Fig. [2ji, 
surprising in the lowest field regime, for the NbB5 sam- 
ple, it is evident the presence of a new cross-over at a 
magnetic field value B cr \ « 9 mT, below which v* de- 
creases with decreasing the external field. In Fig. [2)3, for 
the NbA02 sample, the v*(B) dependence is displayed in 
the low field range for both positive and negative field val- 
ues. In this case it is clear that the critical vortex velocity 
goes to zero crossing the zero magnetic field value, as well 
as the v*(B) curve is symmetric with respect to zero field. 
We outline that all the measured samples show the same 
remarkable feature: the presence of a new cross-over field 
P cr i in the low field region, regardless the thickness of 
the Nb film and the strip width. However, a study of 
this unusual low field behavior on the same Nb film for 
different line widths of the strip is presented in Fig. [3j 
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FIG. 1. I-V curves for different values of B in the low field 
range, at T = 4.2 K, for NbAl sample. The field values 
are: a=0, b=0.78, c=1.56, d=2.08, e=2.34, f=2.60, g=3.90, 
h=5.20, i=6.50, 1=7.80, m=9.10, n=10.40 mT. Inset: Full 
view of the I-V curves where the transition to the normal 
state is shown. 



It is shown that the B cr i cross-over value is shifted 
to lower fields as the line width becomes smaller. This 
means that to be able to detect such experimental be- 
havior, not only a fine variation of the external magnetic 
field should be ensured, but also the line width should be 
sufficiently large. Moreover, in principle, no upper limit 
to the line width can be established. This clearly states 
a size dependence of the cross-over low field. 
The temperature behavior of the B cr \ cross-over field has 
also been investigated. In Fig. H^i the evolution of the 
v*(B) for different temperatures up to T c is reported for 
the NbAl sample. A threshold temperature t* = 0.75 
is found, below which the anomalous low field behavior 
of v* (B) is preserved, with a small variation of the B cr i 
value, between 2.5 and 3 mT. On the contrary, as shown 
in Fig. HJd, above t* the B cr \ crossover field disappears, 
leading to a monotonous decreasing function of v*(B). 
Similar results have been obtained on other samples. 25 
In order to obtain more information on the observed 
anomalous low field behavior v*(B), we visualized the 
first flux penetration by MOI technique with an indica- 
tor filmiii A description of the experimental set-up can 
be found elsewhere^ The results of the measurements 
can be here summarized: 

1) the first penetration of vortices occurs in a discontin- 
uous way, preferably through small circular defects along 
the sample edges; 

2) no dendritic pattern is observed with both transport 
current up to 40 mA and applied magnetic field up to 
80 mT; 

3) the flux penetration is characterized by a typical fan- 
like shaped pattern, as shown in Fig. [5j which becomes 
smoother and smoother as the vortices diffuse deeper in 
the film. 

In addition we have studied the morphology of the Nb 
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FIG. 2. a) Vortex critical velocity v* as function of B in the 
low-field range, at T = 4.2 K for the sample NbB5. The arrow 
marks the crossover field B cr i- b) The v*(B) dependence at 
T — 4.2 K for the sample NbA02 in the low- field range for 
positive and negative field values. 



strips by a FESEM analysis, which confirms the pres- 
ence of circular defects both along the edges and the in- 
ner part of the strips, randomly distributed, whose size 
results in few tenth of micron width. We estimate that 
micro-pores width is around 0.4 /im, and their spacing 
is almost 5 /im, within a cluster of around 20 /im, then 
clusters are randomly distributed in each strip. 



III. DISCUSSION 

The low field dependence of v*(B) can be understood 
taking into account the non uniform magnetic flux pene- 
tration, with the protrusion formation as the field starts 
to penetrate in the mT range (see Fig. [5j). In fact the 
fan-like magnetic flux penetration mirrors the non ho- 
mogeneous state of the fairly strong pinning Nb strips, 
as already outlined in the literature i 16 i 18 
Moreover, by applying the driving current, in the field 
penetrated regions a flux flow state of the moving vor- 
tices is achieved. MOI and direct transport measure- 
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FIG. 3. Critical vortex velocity as function of B in the 
low field region, at T — 4.2 K, for three different strips 
on the 60 nm thick Nb film. Open circles refer to NbAl 
(w — 20 /im) , triangles to NbA3 (w — 50 /xm), and squares to 
NbA4 (w — 100 /im). The arrows mark B cr i for each sample. 



merits give us complementary information: the former 
states the formation of channels due to the fan-like pene- 
tration, the latter states that vortices initially penetrated 
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FIG. 4. Critical vortex velocity as function of B in the low 
field region for NbAl, at different t = T/T c a) At t = 0.5 (full 
circles) and 0.6 (stars), b) At t = 0.90 (open circles) and 
t = 0.99 (full triangles). 




FIG. 5. MOI measurements of the sample NbA4 at T = 4.1 K. 
After zero field cooling, a sequence of increasing magnetic field 
was applied. Dark contrast represents the zero-field Meiss- 
ner state. The bright line is the sample edge along the strip 
length. A fan-like penetration pattern clearly develops after 
vortex nucleation occurred at sample edges. 



in the channels, can be driven afterwards at high veloci- 
ties. Therefore the inhomogeneous magnetic penetration 
induces preferential stationary channels for the moving 
vortices with a total effective width I. In particular, from 
FESEM analysis, these channels can be assumed around 
0.4 fim wide, comparable to the micro-pores, grouped 
together in clusters along the strip. This means that 
the voltages measured across the distance L between the 
voltage tips are actually produced across an effective dis- 
tance I < L; so that in the low field range, reasonably for 
B < B cr i, the flux flow channels fill only partially the 
whole strip length. 

For this reason an analogous scale factor should affect 
all the voltage measurements and in particular the vor- 
tex critical velocity and the flux flow resistance. On the 
contrary when the flux penetration becomes smooth and 
regular Bean-like, above the threshold temperature £*, 
the pinning effectiveness decreases so the vortex flow be- 
comes uniform as well. Therefore at temperatures below 
t* and in the mT range, the flux flow resistance Rff 
is expected to assume lower values than the Bardeen- 
Stephen prediction (RBS), indeed this is shown in Fig.[6k. 
Our hypothesis is sustained by the observed linear de- 
pendence of Rff above the same threshold temperature 
t* ~ 0.75T C . Linear fits of the I-V curves close to V* 
have been performed in order to get the dynamic resis- 
tance in the flux flow state, Rff, which are displayed in 
Fig. |BJl. Inset of Fig. [6^ reports the scale factor L/l as 
function of 5, as deduced by our data. The scale fac- 
tor displays a maximum in the same low field range in 
which the flux flow resistance Rff shows a minimum (see 
Fig. [6^0 . Since Rff turns equal to Rbs — RnB / B C 2 (B C 2 
has been determined with the criterion of J C (B C 2) = 
from the corresponding J C (B)), L/l becomes equal to 1 
around B w 2 mT, which is the same cross-over field B cr \ 
found in the corresponding (B) curves of Fig. [4^. 
In Fig. [6}} we report at t = 0.67 the magnetic field de- 
pendence of v* corrected by the L/l factor for the NbAl 
sample. At low fields, v*(B) is still well different from 
the usually observed power-law dependence v* w B~ x l 2 . 
Nevertheless the very low field dependence of v* has not 
yet been studied in the literature so far. In this field 
range, in the case of isolated vortices, a value of the crit- 
ical velocity for B towards zero, 1$, larger than the LO 
prediction should be found. As the field increases from 
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FIG. 6. a) Flux flow resistance as function of B in the mT 
range for NbAl, at t = 0.67 (squares) and t = 0.74 (trian- 
gles). The corresponding lines refer to the Bardeen-Stephen 
prediction Rbs = R n B/B C 2. Inset: Scale factor L/l as func- 
tion of B at t = 0.67 (solid line) and t = 0.74 (dotted line), b) 
Vortex critical velocity as function of B at t — 0.67 for NbAl 
corrected by the factor L/l (open circles). The full circles 
represent the experimental data. 



zero up to values such that the vortex lattice constant 
ao = (^if) 1//2 is larger than the quasi-particle diffusion 

length l e = (Die) 1 / 2 , the damping coefficient r] of the 
moving vortex lattice is equal to the one for isolated vor- 
tices. For this reason v* should approach an almost con- 
stant value Vq. Only for larger fields, i.e. for ao < l e , as 
for as B < B cr the power-law B~ x l 2 dependence should 
be recovered; finally for B > B cr the LO prediction is 
followed^ In our case we evaluate the quasiparticle dif- 
fusion coefficient D = {\vple) — 2.0 x 10 _4 m 2 /s (using 
for the Nb Fermi velocity the value vp = 2.73 x 10 5 m/s 
[26j and deducing the electronic mean free path, Z e , 
from the measured low temperature resistivity [27j) and 
r e = 5 x 10 -10 s fitting eq. (I), 19 so that the condition 
ao ~ l e corresponds to B* « 39 mT. This means that 
in the low-field range considered in this work neither the 
LO nor the B~ x l 2 predictions should be expected. Nev- 
ertheless, to give a quantitative estimate of the low field 



cut-off vortex critical velocity Vq, the contribution from 
the fan-like penetration should also be added. Indeed 
from MOI measurements it is evident that vortices en- 
tering the sample have a non-vanishing velocity compo- 
nent perpendicular to the Lorentz force (i.e. along the 
direction of the applied bias current), thus leading to 
lower local values of the longitudinal component of the 
flux flow critical velocity v*. Therefore, the almost con- 
stant value found experimentally for B towards zero in 
Fig. [6)3, can be further corrected to the actual value of 
Vq. Moreover, for a quantitative analysis the non uniform 
self field generated by the bias current and local inhomo- 
geneities within the channels should also be considered. 
On one hand, a dramatic influence of self-field effects, 
which could be comparable with the lower applied fields, 
can be excluded. In fact for our samples at low mag- 
netic fields the depinning critical current I c has a steep 
decrease as function of the external B which therefore 
is dominant with respect to the self field. On the other 
hand, a more realistic non constant width of the vor- 
tex channels determines a non-homogeneous distribution 
of vortex velocity. As a consequence locally the critical 
vortex velocity can be reached before the mean critical 
vortex velocity, so this implies a further reduction of the 
measured threshold voltage V * . 

In the end, a possible contribution of vortex- ant ivortex 
(V-AV) pairs to the vortex critical velocity at very low 
fields should be also considered. In particular such an 
effect becomes dominant at B = and close to T c in 
the presence of a current biasing. In principle V-AV 
pairs formation exist, but to be significant a saturation 
of the critical voltage V* to a finite value at low fields 
should occur. By comparing our data on several Nb 
samples, we can reasonably state that no V-AV pairs 
contribution to v*(B) is usually observed, see for exam- 
ple Fig.s [2^ and [2)3, where the critical velocity v* ap- 
proaches zero value crossing the B = axis. Nevertheless 
in Fig.s [4] and [6)3 an almost constant value Vq is found 
for B ~ 0, so that V-AV pairs may contribute at least in 
the case of the smallest line width investigated. 
Finally, we emphasize that, very recently, the observation 
of this new cross-over effect has been confirmed in Al thin 
films grown on top of a periodic array of magnetic micro- 
loops, 28 if magnetized in a flux-closure state.— In this 
case the magnetic state of the artificial pinning structure 
is also effective to sustain a channel-like vortex motion^ 
A comparison with the rather strong intrinsic pinning 
Nb superconductor outlines that pinning, intrinsic or ex- 
trinsic, significantly affect the vortex critical velocity in 
type-II superconducting thin films. 

In summary a major cause of the observed reduced flux 
flow resistance at low applied magnetic fields is very likely 
the inhomogeneous flux penetration due to the micro- 
porosity, which is mirrored in the unusual low field be- 
haviour of the vortex critical velocity v*(B). 
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IV. CONCLUSIONS 



We observe a low field cross-over effect in the mag- 
netic field dependence of the critical vortex velocity in 
wide Nb strips. We believe that this cross-over effect in 
the low- field regime can be justified by vortex channeling 
through the inhomogeneous pinning center distribution, 
as emphasized by MOI measurements. We observe that 
B cr \ increases by both increasing the line width and the 
thickness of the samples. This suggests that such cross- 
over effect can be hidden by the sample size. 
The channeling effect on the moving vortices is fur- 
ther confirmed by the magnetic field dependence of the 
flux flow resistance which turns back into the Bardeen- 
Stephen linear behavior as the B cr i field is overcame. 



Therefore B cr \ marks the crossover between the channel- 
like vortex motion and the uniform flux flow state. This 
observation is likely to be a common feature in type-II 
superconducting thin films exhibiting moderately strong 
pinning, as preliminary observed in the case of Al thin 
film grown on artificial pinning structures^ and here 
clearly established for the intrinsic Nb superconductor. 
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